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Abstract 

We report on parametric amplification in dynamic radiation force produced by a bichromatic 
acoustic beam in a fluid. To explain this effect we develop a theory taking into account the 
nonlinearity of the fluid. The theory is validated through an experiment to measure the dynamic 
radiation force on an acrylic sphere. Results exhibit an amplification of 66 dB in water and 80 dB 
in alcohol as the difference of the frequencies is increased from 10 Hz to 240 kHz. 
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Acoustic radiation force in fluids is a phenomenon that has been investigated for over 

i 

a century It results from variations in energy and momentum of the wave as a conse- 
quence of scattering, attenuation, or distortion by nonlinear effects. This force is similar to 
the optical radiation pressure exerted by electromagnetic waves on responsive particles |2|. 
Radiation force can be either static or dynamic with respect to its time dependency. Static 
radiation force is a time-averaged quantity produced by a monochromatic wave that corre- 
sponds to the dc component in the spectrum of the stress (pressure) [3| . Dynamic radiation 
force can be generated by a bichromatic acoustic wave. This force is associated to the 
difference frequency component in the spectrum of the stress 

One of the first applications of dynamic radiation force was devised in 1928 by Sivian \^ 
to measure acoustic power on a suspended disk. In 1953, Macnamara et al. [Of developed a 
method based on this force to measure absorption in liquids. After these pioneer applications, 
dynamic radiation force has passed unnoticed until late 1970's. In the last three decades, it 

Q 

has been applied for measuring ultrasound power of transducers 0], inducing oscillation in 
bubbles |J and liquid drops Q, and exciting modes in capillary bridges [l^]. Furthermore, 
dynamic radiation force is the underlying principle of some elastography imaging techniques 



such as shear wave elasticity imaging 



llj and vibro-acoustography 



Despite early applications, dynamic radiation force was only recently investigated in 
theoretical grounds. Mitri et al. Il3j calculated it on elastic cylinders. He also studied the 

n n 

force produced by a bichromatic standing plane wave [14]. Silva et al. [4j obtained the 
dynamic radiation force on elastic spheres which was confirmed in an experiment realized 
by Chen et al. [15]. We emphasize these results are only valid when the difference between 
the fundamental frequencies of the wave is very narrow. So far, no one has taken into 
account the influence of the fluid nonlinearity in dynamic radiation force thoroughly. The 
nonlinearity of the fluid is described by the thermodynamic relation p oc p( 1+B / A ) , where p 
and p are respectively the pressure and the density of the fluid, and B/A > is commonly 
used in acoustics as the nonlinear parameter It is worthy to note that static radiation 
force in a ideal fluid does not depend on the nonlinearity parameter B/A \x\. On the other 
hand, one may inquire on how does dynamic radiation force depend on the nonlinearity of 
the fluid? 

In this letter we undertake this question. Our analysis unfolds that dynamic radiation 
force may achieve a regime of parametric amplification. The concept of parametric amplifi- 
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cation arose in radio engineering and is widely known in optics In acoustics, parametric 
amplification can also be understood as follows. The mixture of two waves of differing an- 
gular frequencies uj\ and u 2 {uj 2 > ^i) generates two new waves; one of which has frequ ency 
equal to u)\ + tu 2 , while the other arises with the difference frequency uj 2 i = lu 2 —lui [19, 20]. 
To demonstrate the parametric regime in dynamic radiation force, we calculate this force 
on a rigid sphere taking into account the nonlinearity of the fluid. An experiment using a 
laser vibrometer is designed to measure the dynamic radiation force on an acrylic sphere 
immersed in degassed water and ethyl alcohol. Results show an amplification of this force 
of up to 66 dB in water and 80 dB in alcohol as the difference frequency varies from 10 Hz 
to 240 kHz. 

Consider a homogeneous and isotropic fluid with adiabatic speed of sound Cq, in which 
thermal conductivity and viscosity are neglected. The fluid has infinite extent and is charac- 
terized by the following acoustic fields: pressure p, density p, and particle velocity v = — V0. 
The function is the velocity potential and V is the gradient operator. These fields are 
function of the position vector r and time t. At rest, these quantities assume constant val- 
ues p = po, p = po, and v = 0. The acoustic fields are governed by the dynamic equations 
of ideal fluids. By using the regular perturbation technique, one can expand the velocity 
potential in terms of the Mach number e <C 1 as <p — £<fi^ + s 2 4>^ + 0(e 3 ). Any analysis of 
radiation force has to be done considering at least the second-order terms of this expansion. 

The excess pressure in the fluid can be written as p — po = p^ +p^ + 0(e 3 ), where j^j 

(i) m 

and 



2c 2 \ dt J 2 dt 



are the acoustic and the second-order pressure fields, respectively. The first two terms in 
the right-hand side of Eq. (j2j) depend only on They correspond to the Lagrangian 

density of the wave. As we shall see the potential does not depend on the nonlinearity 
of the fluid, while (p^ does. Thus, Eq. (J2J) has contribution of two terms called here the 
Lagrangian and the nonlinear pressures. 

Let So be the surface of the object target at rest. One can show that the instantaneous 
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force on the object up to second-order in the excess of pressure is given by 



p {1) ndS- // p {2) ndS, (3) 

'S(t) J J S 

where S(t) is the moving object surface, n is the outward normal unit-vector on the integra- 
tion surface. Assuming that the sphere is under influence of a bichromatic acoustic beam 
with fundamental angular frequencies uj\ and U2, the dynamic radiation force is produced 
by the contribution of stresses at the difference frequency uj 2 i — uj 2 — 

Consider the Fourier transform of a function g(t) as and its inverse denoted by 
. The dynamic radiation force is given in terms of the component of Eq. (jHJ) at uj 2 \. 



Accordingly, the dynamic radiation force is |4| f 2 i = T 



--1 



121 



where 



1*21 = / / n • S 21 dS - iepQ<jj 2 \T / / (fi^ndS 
J J s JJs(t) 



where i is the imaginary unit and 

S 21 = -^[p( 2 )| + p oV «v«] 



(4) 



(5) 



is the amplitude of the dynamic radiation stress with I being the unit tensor. The dyad 
Pov^-M 1 ) is the Reynolds' stress tensor. 

To obtain the dynamic radiation force over an object, we have to solve the corresponding 
scattering problem described by and (f)( 2 \ These functions satisfy the linear and the 



second-order wave equations 22 

□ 2 0« 

□ 2 0( 2 ) 



]_d_ 
c 2 dt 



(6) 
(7) 



where D 2 = V 2 — (1 / 'cq)(<9 / 'dt) 2 is the d'Alembertian operator and 7 = 1 + B/A. Hence, the 
scattering problem should be solved through Eqs. (0) and (JJJ) with appropriate boundary 
conditions. 

Now consider a bichromatic plane wave formed by the excitation in the velocity field 
v = eco{sm.uj\t + sinu^)^ at z — 0, where e z is the unit vector in the z direction. For the 
perturbation terms of the velocity potential, one has the boundary conditions d<fP~> /dz = 
— Co(sinu;it + sinco> 2 t) and d<p^ jdz = both at z — 0. The solution of Eq. © is 



c Re \ L e -^t-kiz) + }_ e -i{u2t-k*z) 
h k 2 



(8) 



where Re means the real-part of a complex variable, ki = uj\/cq, and k 2 = oj 2 /cq. We are 
only interested in the second-order velocity potential at the difference frequency u 2 x ■ Hence, 
from Eq. (J7J) we have 



where /c 2 i = ^21 /co- 
in what follows we calculate the dynamic radiation force on a rigid sphere of radius a. The 
amplitude of the force in Eq. (J1J) has contributions from the Lagrangian, the acoustic and the 
nonlinear pressures, and the Reynolds' stress tensor. One can show that the contribution of 
the second term in the right-hand side of Eq. (j3J) is proportional to UjmZ m )~ l , m = 1,2. The 
quantity Z m is the mechanical impedance of the oscillating sphere 15[ at the frequencies uo\ 
and uj 2 . In the experimental setup we use frequencies above 2.2 MHz; hence, the magnitude 
of f small as 10 6 , which is much smaller than other contributions to the 

dynamic radiation force. Therefore, we neglect this contribution here. 

Let us focus on the scattering of the nonlinear pressure at the difference frequency p 2 i — 
e 2 pod^i 1 'dt. We restrict our analysis to the case k 2 \z ^> 1 and a Zq. The energy density 
at the acoustic source is Eq = e 2 poc 2 ] /2. From Eq. @, the amplitude of the nonlinear 
pressure in the spherical coordinates (r, 9, ip), is given in terms of partial spherical waves as 

+00 

pa! = A^(2n + l)i n \j n {k 21 r) + b n h n {k 21 r)]P n {cos9), (10) 

n=0 

where A = —iE 'jk2iZoe' lk21Z0 , the functions j n and h n are the spherical Bessel and first-type 
Hankel functions of nth-order, respectively. The function P n is the Legendre polynomial of 
nth-order. The scattering coefficients are given by j3] 

b n = -j' n (x 2 i)/h' n (x 2 i), n^l 

(P0/Pl)jl(«2l) - X 21 j[(x 2l ) 



(po/pi) h(x 21 ) - x 2 ih' 1 (x 21 ) 
where x 2 i = k 2i a. The symbol ' stands for the derivative. 



By substituting Eq. (jlO)) into Eq. 
pressure, and using the result of Ref. 



neglecting the direct contribution of the acoustic 
one obtains the dynamic radiation force as 



f 2 i = na 2 E Y 21 e- i ^ t - k ^ zo) e z . (11) 
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The dynamic radiation force function Y 2 \ is given by 



Y21 = - — {—{R* 11 R2,2 + Rl 2 R2,i)+i2(n + i) 
xix 2 L pi • ~o 

x [xix 2 - n(n + 2)] (R\ n R^ n +i + Rl, n +i R 2,n) J 



(po/Pl)^l(^2l) - ^21^i(^2l) ' 



where 



Rm,n 2 hi I \ ' ^ / 1) 



a^m = with m = 1, 2. The symbol * is the complex conjugate. The first term of Eq. (j!2|) 



is the dynamic radiation force used in the literature |4J, |13j]. The last term in the right- 
hand side of Eq. (|12|) is due to the nonlinearity of the fluid and is associated to parametric 
amplification in dynamic radiation force. This term has not been treated in previous works. 
The regime of parametric amplification can only be neglected when io\ ~ cu 2 . It is important 
to note this amplification depends on the quantities, namely, the ratio zo/a, 7, and ^21- 

To verify Eq. (fTT|) . we realized an experiment to measure dynamic radiation force. The 
basic experimental setup is described in detail in Ref. In Fig. [T]we illustrate the ex- 
perimental apparatus. An in-house flat transducer with diameter d = 22 mm is used to 
insonify an acrylic sphere inside a tank filled either ethyl alcohol or degassed water. The 
sphere has radius a = 1.6 mm, density pi = 1190 Kg/m 3 and is suspended in a pendu- 
lum along the beam axis 20 cm away from the transducer. The parameters for alcohol 
(water) are p = 785 (1000) kg/m 3 , c = 1100 (1500) m/s, and 7 = 11 (6). We used a 
suppressed-carrier amplitude-modulated signal to drive the transducer. The carrier fre- 
quency is Uo/2tt = 2.25 MHz and the modulation frequency c<j 21 /27r is swept from 10 Hz to 
240 kHz. In such configuration the wave frequencies are 00^2 = ujq =F ^21/2. For the specified 
difference frequency range, we measured the acoustic pressure in water at 200 mm away 
from the transducer aligned with the beam axis. The measurement was performed by a 
polyvinylidene fluoride membrane hydrophone (Y-33-7611, GEC-Marconi, Great Britain). 
The amplitude of the measured pressure remained constant within an error of less than 5% 
as the difference frequency varied in the specified range. The dynamic radiation force on 




FIG. 1: (color online) Experimental apparatus utilized to measure the dynamic radiation force on 
the sphere. The dimensions of the tank is 100 x 64 x 38 cm for water and 40 x 40 x 24 cm for 
alcohol. 

the sphere is obtained by measuring vibration velocity of the sphere. The force is given by 
the product of the vibration velocity and the mechanical impedance of the sphere at lo 2 \. 
The vibration of the sphere was detected by a laser vibrometer (Polytec GmbH, Waldbronn, 
Germany), which was aligned with the beam axis of the transducer. The signal from the 
vibrometer was filtered by a lock-in amplifier (Perkin Elmer 7265, Oak Ridge, TN) at the 
difference frequency uj 2 i. 

The frequencies used in the experiment are such that k i2 a ^> 1. Thus the dis- 



tance at which the farfield of the transducer begins for the frequencies lo\ and uo is \24\ 
Zi,2 = &i,2^ 2 /87T, which corresponds to Z\ ~ z 2 = 18 cm. Furthermore, in this region the 
incident wave resembles a plane wave with circular cross-section whose 3 dB-radius is ap- 
proximately [2J] 0.35d/2 = 3.9 mm. We can reasonably assure that the sphere is thoroughly 
inside this region and the plane wave approximation can be used for the incident beam. 

The dynamic radiation force as described in Eq. (fTTJ) is only valid within the preshock 
wave range. When a finite-amplitude wave propagates in a fluid, its form tends to steepen up 
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TABLE I: Summary of the experimental measurements. 



Mach 


Dynamic Radiation Force Gain 


Number 


10 Hz 


240 kHz 


alcohol 1.7 x 1(T 5 


3.0 fiS 


34 mN 80 dB 


water 1.0 x 1CT 5 


2.8 /iN 


6.7 mN 66 dB 



to develop shocks. The formation of shock waves implies dissipation, which is not described 
by Eqs. © and (JZj). For a monochromatic plane wave the preshock range has length 2^ 
i = A[27re(l + B /2A)]~ 1 , where A is the wavelength. The only unknown parameter necessary 
to determine the preshock wave range is the Mach number. However, we may estimate it 
by equating the dynamic radiation force measured at u;2i/27r = 10 Hz to Eq. (JTTJ). In this 
frequency range the theory of dynamic radiation force was already validated experimen- 
tally 15[. A summary of the measured data in the radiation force experiment is shown in 
Table HI According to the measured data, the preshock wave range is I = 62 (272) cm for 
alcohol (water). Therefore, the sphere at z = 20cm is inside this range for both alcohol 
and water experiments. 

In Fig. |21 we exhibit a comparison between the measured dynamic radiation force and 
Eq. (fTT]) . Prior to the vertical dotted line at 2 kHz, parametric amplification has a minor role 
in the dynamic radiation force. This region coincides with the result obtained in Ref. j^|. 
Beyond this line we have a prominent amplification of the dynamic radiation force. The 
theory predicts a gain of 52 (44) dB in alcohol (water) as the difference frequency varies 
from 10 Hz to 240 kHz. Discrepancies between the theory and experimental results might be 
related to diffraction, thermoviscous effects of the fluid, and surface and elastic properties of 
the sphere. Factually, the acrylic sphere is a viscoelastic material which allows the formation 
of internal and surface waves. These waves are not present in a rigid sphere. Thermoviscous 
effects may produce streaming in boundary layers surrounding the sphere and nearby the 
walls of the tank. Streaming can cause a blueshift in the waves near to the surface of the 
sphere; thus increasing the dynamic radiation force. A theory with thermoviscous effects 
requires a full solution of the scattering problem using the Navier-Stokes equations |26j |. 

In summary, we showed that dynamic radiation force is subjected to parametric amplifi- 
cation. Measured amplification were 80 dB in alcohol and 66 dB in water as the difference 
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FIG. 2: (color online) Comparison of the theory and experiment to measured dynamic radiation 
force on the acrylic sphere suspended in water (o) and ethyl alcohol (o). 

frequency was tuned from 10 Hz to 240 kHz. Results are in reasonable agreement with the 
theory presented in this letter. We believe parametric amplification may set new applications 
of radiation force in elastography and trapping particles. 
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